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Abstract. The problem of damping a system of linear oscillators is consid¬ 
ered. The problem is solved by using a control in the form of dry friction. The 
motion of the system under the control is governed by a system of differen¬ 
tial equations with discontinuous right-hand side. A uniqueness and continuity 
theorem is proved for the phase flow of this system. Thus, the control in the 
form of generalized dry friction defines the motion of the system of oscillators 
uniquely. 
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Introduction 

The present work is closely related to the talk of the one of the authors (A.I.O.) 
at the international seminar “Control Theory and Theory of Generalized Solutions 
of Hamilton-Jacobi Equations” (CGS’2015) dedicated to the memory of A.I. Sub¬ 
botin. The subject of this work correlates with one of the central topics of work 
of Andrei Izmailovich: How to define the solution of a problem so that the whole 
corresponding theory takes an attractive and ultimate form. 

Of course, Andrei Izmailovich and we deal with substantially different problems. 
Our work is devoted not to solution of nonlinear partial differential equations but 
to ordinary differential equations (ODEs) and related linear transport equations. 
Nevertheless, the basic idea of the work consists in obtaining the final existence and 
uniqueness of the motion for a very specific system of differential equations with 
discontinuous right-hand side (RHS) on the basis of a nonconventional concept of 
a solution. 

In the case under consideration, the system arises in an attempt to control in 
a quasioptimal way a system of an arbitrary number of linear oscillators by using 
a feedback control in the form of a generalized dry friction [1-3], and the adopted 
concept of solution is suggested by the classical work by R. DiPerna and P. Lions [4] . 

1. Problem statement and background 

As it is well-known, by using the maximum principle one can explicitly construct 
a feedback control for the minimum-time damping of a single linear oscillator [5]. 
A direct generalization of such task is the problem of damping a system consisting 
of an arbitrary number N of linear oscillators with different eigenfrequencies oJi: 

(1.1) X = Ax + Bu, xGY = M.'^^, u S U = R, juj < 1, 
where the matrix A and the vector B have the form 

(1.2) A, = J), A = diag(A,), B = ®B,. 
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In the natural coordinates {xi,yi), i = the system can be written as 

follows: 

Xi = Hi, 

yt =-ujfxi + u, |m| < 1 , i = 

Asymptotical properties of the system are determined primarily by the presence 
or absence of resonances, i.e., non-trivial relations between eigenfrequencies of the 
form 

N 

(1.3) = 0, 0 7 ^ m = (toi, ...,toat) G 

i=l 

The Kalman controllability criterion [ 6 ] for system (1.1)-(1.2) reduces to the dis¬ 
tinction of the eigenfrequencies: uji ^ ujj at i ^ j and it is, of course, substantially 
weaker than the condition of absence of resonances. 

The minimum-time problem for system (1.1)-(1.2) can be reduced to the boundary- 
value problem of the Pontryagin maximum principle 

X = Ax + Bu, p = —A*p, 

(1.4) 

u = sign(i?,p), x(0) = xo, x(T) = 0, h(x,p)=0, 

with respect to the Hamiltonian 

h{x,p) = {Ax,p) + \{B,p)\ — 1 = max{(Aa:,p) -I- {Bu,p) — 1}, 

|u|<l 

where angle brackets (.,.) stand for the standard scalar multiplication in | • | is 
the Euclidean norm, and the maximum is taken over the interval {u G R : |u| < 1}. 

We note that system (1.4) is the Hamiltonian one with 2N degrees of freedom 
and -|- 1 integrals of motion. These integrals are the Hamiltonian h and energies 

li = \ ijii + i = 

of normal oscillations of the vector p written in the from pi = where the 

variables and pi are dual to Xi and yt, respectively. From the point of view of 
the canonical system of the maximum principle (1.4), the problem of the optimal 
damping of a single linear oscillator is completely integrable, because in this case the 
number iV -|-1 of integrals of motions is equal to the number 2N of degrees of free¬ 
dom. Note that the similar identity is the basic assumption of the Liouville-Arnold 
theorem on the complete integrability of Hamiltonian systems [7]. In contrast, prob¬ 
lem ( 1 . 1 )-( 1 . 2 ) is probably not completely integrable, so that an analytic design of 
the optimal control by using methods based on the Pontryagin maximum principle 
is hardly possible. 

The general problem to be solved is to construct a nonoptimal feedback control 
steering the system to equilibrium. An approach to design of an asymptotically 
optimal feedback control for system (1.1)-(1.2) was suggested in [1,2]. When using 
methods of [ 1 , 2 ], the ratio of the motion time to zero under the suggested control 
to the minimum time is close to one provided that the initial energy of the system 
is sufficiently large. 

Within the framework of the present work, we confine ourselves with issues 
related to a construction of the control in the generalized dry friction form and 
to study of the motion of the system under this control. We describe an approach 
to design of the control, study systems of differential equations with a discontinuous 
right-hand side related to the motion of the system, and show that the motion can 
be defined uniquely. In terms of Refs. [1,2], this gives a description of the dynamics 
of the system under the asymptotically optimal control within the high-energy 
regions of the phase space. In the present work, we provide proofs of existence 
and uniqueness theorems for arising deferential equations with discontinuous RHS. 
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Investigations of asymptotic properties and detailed description of the control can 
be found in [2,3] 

It should be noted that substantially different approaches to construction of 
the feedback control for linear systems are known as well, for example, the methods 
based on the Kalman approach to the feedforward control [8,9]. Papers [8,9] include 
estimates for the motion time under this control. The time is comparable with the 
optimal one: the ratio of the motion time to zero under this control to the minimum 
one is bounded. 


2. Control for a system of oscillators 

A well-known geometric interpretation of the Hamilton-Jacobi-Bellman equa¬ 
tion for the minimum-time problem is that the momentum vector dT/dx at a point 
X is an inner normal to the reachable set T>[T(x)), where T{x) is the optimal time 
for the controllable system. 

Definition 1. A reachable set 'D{T) is the set of ends of admissible trajectories of 
a controllable system from zero and parametrized by the time interval [O.T]. 

Optimal control has the form 

dT 

(2.1) u{x) = -sign{B,p{x)), p = 

where p is the outer normal to the reachable set 'D{T{x)) form zero, whose boundary 
passes through x. 

The control in the dry friction form 

N 

(2.2) M =-sign^ 

i=l 

where \i are some positive coefficients, depending on positions and velocities, 
arises when an asymptotic approximation is substituted for the exact reachable 
set 'D{T{x)) of system (l.l)-(1.2). According to the asymptotic theory of reach¬ 
able sets [10-12], as time T —>■ oo TH, where H is a fixed convex body, serves as a 
good approximation to X>(T). The body H can be uniquely defined by its support 
function. 


Definition 2. The support function Hm{0 of ^ closed convex set M has the form 

(2.3) HMiO = svcp{^,x) 

x^M 

and defines the set M uniquely [13]. 


More precisely, for the system (1.1)-(1.2) under consideration. 


Theorem 1. Suppose that a momentum p is written in the form p = 
where is the dual variable for Xi, pi is the dual variable for yi, and Zi = [pf + 
In the nonresonant case, when there are no nontrivial relations between 
eigenfrequencies of the form (1.3), the support function Ht of the reachable set 
'D(T) has as T ^ oo the asymptotic form 


(2.4) 


Ht{p) = T 



Zi cos (fii 


dp + o{T) = TS^{z) -I- o(r), 


and the support function of the convex body H is given by the main term f){z). 
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Proof. By definition, the support function of the set 2?(T) has the form (2.3), where 
the supremum is taken over admissible controls and x{T) is the state at time T of 
the control system (1.1)-(1.2) such that a;(0) = 0. 

By using the Cauchy formula, we conclude that 


fT pT 

{x{T),p) = / Bu{t),p)dt = / u{t)B*e^ ^'^~*^pdt. 

Jo Jo 


Taking the supremum under the integral sign and the change of variables t i-A T — t, 
we obtain 


(2.5) 



fT 

sup u{t)B*e^ ^'^-*'>pdt = / 
|ii(t)|<l Jo 


\B*e^‘^p\dt. 


In the dual coordinates formula (2.5) takes the form 


H-D{T)ip) — 

Jo 


N 


rji cos Wit + Wj sin Wit 


dt. 


We note that this expression represents an integral of the function 




N 

Vi cos (fi + sin ip^ 


taken over the rectilinear winding (pi(t) = ojit of the torus T = (K/27rZ)^ with 
angular coordinates ipi. 

Suppose that there are no resonances in the system, i.e., there are no nontrivial 
relation between eigenfrequencies of the form (1.3). Then the time average coincides 
with the space average [7]: 

lim ^ f f{ip{t))dt= f f{p)dip 

T-^oo 1 Jq Jj- 

We note that 

r]i cos ipi + w“^^iSin(/Ji = ZiCOs{ipi + at), 
where a = (a^) is a constant point of the torus. Therefore, 


[ f{p)dp= f f{p-a)d(f= f 

j 7” 77 ” 77 ” 

Thus, we obtain the statement of the theorem 

(2.6) _^lin^ ^HviT){p) = 


N 


Zi cos Pi 


dp. 


IT 


N 

Zi cos Pi 


dp. 


□ 


Note that the theorem follows from the general theory describing the asymptotic 
behaviour of support function of reachable sets of linear systems developed in detail 
in [11]; however the above proof is much simpler than the general theory. 

The support function Hq{p) of the convex body Q is the main term of the 
asymptotic expression in (2.4) 


Hn{p) = ^{z) 



Zi cos Pi 


dp, where z 


[zi,..., zn) € 
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The vector p is normal to the boundary of dft at the point dHQ{p)/dp. The normal 
vector to the approximate reachable set pil, whose boundary passes through x, is 
defined by the equation 


(2.7) 


_ dHa{p) _ dS^[z) dz 
dp dz dp ’ 


where p € and p > 0 are unknown. The function is differentiable, and Eq. 
(2.7) has unique solution, because the boundary of n is smooth [12]. The strategy 
of the control design, which follows from Eq. (2.7), can be applied in the resonant 
case as well, when the asymptotic expression (2.4) does not work; however in this 
case quasioptimal properties of the control are lost. The function p = p{x) from 
Eq. (2.7) plays the same role for the considered control as the optimal time T{x) 
for the optimal control (2.1). The momentum p in (2.7) has the form p = dp/dx. 
The function p{x) is the norm of the vector x in the metric such that the body 17 
is the unit ball. This is a smooth function outside zero. 

Although using the control in the dry friction form (2.2) helps to damp oscilla¬ 
tions, it does not necessarily lead to a complete stop of the system. Strictly speaking, 
standstill zones might appear, where the system is not moving at all, despite the 
fact that the equilibrium is not yet reached. The method suggested in [1-3] com¬ 
bines several strategies of control applied successively at high, medium, and low 
energies. At high and intermediate energies, the scalar control in the form of the 
generalized dry friction is used (2.2). At low energies, a substantially different feed¬ 
back control law, which is constructed by using common Lyapunov functions [1], 
is employed. In the present work, we are interested in the motion within regions of 
high and medium energies. 


3. Motion under the control 


Differential equation with discontinuous RHS occur naturally in optimal control 
theory. A conventional approach to the problem of existence of solutions for such 
equations is based on the Eilippov theory of differential inclusions [14]. However, 
an intuitive concept of a controlled motion includes both an existence and a unique 
determination of the trajectories of the system by a control law. The corresponding 
uniqueness issue is generally beyond the Filippov theory. 

The control in the generalized dry friction form also leads to the motion of a sys¬ 
tem. This motion is formally described by a differential equation with discontinuous 
RHS: 


(3.1) 


X = Ax — B sign 



u{x) = —sign 



Here u{x) is a multivalued function, because sign(O) is defined non-uniquely and can 
take any value in the interval [—1,1]. In fact, we deal with the differential inclusion. 

As it turns out, a motion of the system under the control in the generalized dry 
friction form, which is described by the differential inclusion (3.1), can be defined 
uniquely. Toward this end one can apply the DiPerna-Lions theory of singular 
DDEs [4]. 

3.1. DiPerna—Lions theory. If b{x) is a Lipschitz function, then the Cauchy 
problem for the ODE 


(3.2) X = h{x), x(0) = xq 

and for the partial differential equation (transport equation) 

(3.3) 


dv h f \ 

= 2 ^ b,{x) 


dt 


dxi'' 


v{x, 0) = vo(x) 
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are equivalent. The method of characteristics says that the solution v of problem 

(3.3) is given by the formula 

(3.4) v{x,t) = VQ{(t)t{x)), 
where (pt is the phase flow for (3.2). 

In the paper of DiPerna and Lions [4], the Lipschitz condition dhjdx € Loo 
is substantially relaxed. Instead of it, the Lipschitz condition in an integral form 
db/dx € Li is imposed. It is shown that the solution of problem (3.3) still exists, 
is unique, and is given by formula (3.4). Thereby, is was demonstrated that one 
can efficiently work with a differential equation whose RHS satisfies the Lipschitz 
condition in the integral sense rather than pointwise. 

The DiPerna-Lions theory is based on the notion of renormalized solution. 

Definition 3. A weak bounded solution v of the Cauchy problem (3.3) is called 
renormalized solution, if for any smooth function /3 : R —>■ R function /3(u) is again 
a weak solution. 

3.2. Motion under the generalized dry friction control. The idea behind the 
DiPerna-Lions theory is to construct a global phase flow, perhaps not everywhere 
uniquely defined, instead of a solution of an individual Cauchy problem for every 
initial condition. Our main result claims that in the phase space of system (3.1) we 
can define a semiflow, which is continuous, uniquely defined everywhere, and gives 
a solution of the form (3.4) to the transport equation. 

Theorem 2. There exists a continuous semiflow x i-A x{t) = 4>t{x), < > 0 such 
that v{x,t) = v{4)t(x)) is the unique renormalized solution of the Cauchy problem 
for the transport equation 

(3.5) ^ = (^Ax - Bsign(^B,^{x)^ , v{x,0)=v{x). 

Each curve t i—>■ x{t) is absolutely continuous, and the differential inclusion (3.1) 
holds. 


Proof. We confine ourselves to a proof of existence of a continuous bounded solution 
of the transport equation (3.5), which is obtained as a limit of classical solutions of 
regularized equations. The remaining statements can be proved by using standard 
techniques from [4,15]. 

We use a two-parameter approximation of the problem. First, we choose a pa¬ 
rameter n —>■ oo such that smooth convex functions m„ : R —>■ R uniformly ap¬ 
proximate the function x i—>■ |a;|. Then, the derivative s„ = rn(, approximates the 
function sign(x) in Li. Note that xsn(x) > 0 for any a; € R. Second, we choose 
another parameter <5 0, so as to freeze the motion under system (3.1) within the 

i5-neighborhood Us = {p{x) < 5} of zero with respect to the distance p. In other 
words, we approximate ODE (3.1) by the nonsingular equation 

(3.6) X = Ax - Bsni^B,^^^^ 


in the domain Vs = {x € R^^ : p{x) > (5}. It is important that all the neighbor¬ 
hoods Us are invariant under the phase flow of (3.6) for positive times because the 
radius-function p is nonincreasing along the phase trajectories. Indeed, the following 
inequality holds: 




< 0 . 


P — 
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There is a duality relation between the support function H and the function p 
(for details, see [2,3]): 


(3.7) 


d'^H d^p dp dH 


dp'^ dx^ dx dp ' 

By using relation (3.7), we rewrite equation (3.6) in the gradient form: 


13 Sj] 


B 


dp 

dx 


= pa{x)-^m„ 


B, 


dp 

dx 


+ XSn 


B 


dx 


B 


The latter can be regarded as an approximation to 




B 


’ dx 


+x 


B 


’ dx 


a{x) = 


d'^H 

dp^ 


dx 


H = Ha. 


In particular, the ODE takes the form: 


0 , 


(3.8) X -- 

The functions involved 
f{x) = Ax — xs 


Fix) = fix) - gix)-^mn iHx)) , if a; is in 14, 


B 


dp 


B 


’ dx 


if X is in Ug. 


g = pa, h = (B 


dp 

dx 


are rather smooth: they are locally Lipschitz outside zero. Eqs. (3.2) approximates 
Eq. (3.1) rewritten in the gradient form 


X = Fix) = fix) - gix)-^\hix)\, fix) = Ax - x 


B 


dp 

dx 


It is important for us that the matrix g = pa is symmetric and nonnegative. By 
omitting the subscript n, we find that the corresponding transport equation takes 
the form 

dv 

= fiVi - gijhjVtsih) = FiVi, 

where 

d d 

^ ^ ^ signli, Fi = fi- gijhjsih), 

and the Einstein summation notation is used. By differentiation, we obtain the 
following equations for the vector-function V with components Vk- 

dv 

(3.9) — = FiVk,i + fi,kVi - gij,khiVisih) - gijhjkVisih) - g^jhjhkViSih), 

where Vk,i = |^, hjk = gij,k = and 5 = 5n denotes m". Equation (3.9) 

is again a transport equation with extra terms fi,kVi — gij^khiVisih) — gijhikVisih) — 
gijhihkViSih) in the RHS. Fortunately, the most “dangerous” and singular term 
Cfc = gijhihkViSih) has a positivity property: 


gkivicTk = gklhkVlg^jhjVi5ih) = (^gkihkVi^ Sih) 


is a positive measure. 


All the other terms are linear functions of V with coefficients bounded outside 
any neighborhood of zero. This implies that w = {gV, V) = gkiViVk is a kind of 
quadratic Lyapunov function: 

(3.10) ^<F,yj^ + LW, w=\v\^ = Y,yl 

Here L is a function uniformly bounded outside any neighborhood of zero. Since the 
matrix g = pa is not strictly positive definite, W cannot be estimated via w, and 











8 


ALEXANDER OVSEEVICH AND ALEKSEY FEDOROV 


Eq. (3.10) is insufficient for establishing an a priori bound for w, not to mention 
W. Nonetheless, we can use the estimate 

(3.11) W = Y,vl<c(^i^XkVk)\{gV,V)y 


where C is a positive function bounded outside any neighborhood of zero. The 
bound holds because the kernel of the matrix g{x) is the one-dimensional subspace 
of the phase space, generated by x. In view of inequality (3.11), we have to find an 
estimate for z = '^XkVk = Ev, where E is the Euler operator Ev = By 

applying the Euler operator to the transport equation (3.9), we obtain: 

3z 

(3.12) — = E^Evi + {EEi)vi = FiZi - FiVi + [EFi)vi. 


Here we use the commutation relation 


Aij = £; A + A 

dxi dxi dxi 


which implies that Evi = Zi — Vi. It is easy to compute EFi\ The function 


F{x) = Ax — Bs 



is clearly the sum of the homogeneous functions Ax and —Bs{B, dpjdx) of degrees 
1 and 0. Therefore, EFi is a locally bounded function. Relation (3.12) now implies 
that 

(3.13) ^<F,y, + C'W, 

where y = z^, and C' is a locally bounded function. Inequality (3.11) says that 
bE < C{y + w). Therefore, by summing inequalities (3.10) and (3.13) we obtain 
that 

A < + MY, Y = w + y, 

where the function M is locally bounded outside zero uniformly wrt the scale n. 

Inequality (3.10) is the crucial estimate that enables us to show that the flow 
X i-A- $t(x) = ^n,t{x) corresponding to Eq. (3.2) is locally Lipschitz, and besides the 
corresponding Lipschitz constant does not depend on the parameter n. Therefore, 
by passing to the limit n —>■ oo we conclude that there exists a Lipschitz limit of 
which defines the semiflow (j)t (x) of Theorem 2 within Vs. Since 5 is arbitrary, this 
proves in particular that the map x i—>■ (l>tix) is continuous if x ^ 0 and (fyix) ^ 0. 

It is in fact obvious that the map x i—>■ iptix) is continuous at zero, because the 
flow (j) maps any neighborhood Us of zero into itself. It remains to consider the case 
X ^ 0, (j)t{x) = 0. Put r = inf{t > 0 : (/>t(a^) = 0}. It suffices to show that (jiriu) is 
close to 4)t{x) = 0 if y is sufficiently close to x. We know already that for any e > 0 
the point 4>r-e(x) depends on x continuously. On the other hand, it is obvious that 
the map 1 1 —> (f’tiy) is uniformly Lipschitz for y in a neighborhood of x. Therefore, 


(3.14) |(/)^(y) - (l)r{x)\ < Clel -b |(()r-e(y) - (t)r-e{x)\. 

Since e is arbitrary and |0T-e(y) — 4>T-eix)\ is arbitrarily small if y is sufficiently 
close to X, the continuity is proved. □ 


A similar phenomenon was discovered by LA. Bogaevskii [16] for the gradient 
differential equations x = —df/dx, where / is a nonsmooth convex function. 
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Conclusion 

In our work, the control for damping a system of oscillators in the generalized 
dry friction form was studied. As it customarily happens in optimal control the¬ 
ory, this control gives rise to differential equations with discontinuous RHS. In the 
present work, we demonstrated that for the case under consideration the problem 
of existence and uniqueness of the motion under the control can be resolved in 
the framework of the DiPerna-Lions theory. It seems important to study similar 
problems for the optimal control. 

An interesting development of the considered problem of damping a system of 
oscillators is given by passage to the infinite-dimensional case. For example, the 
problem of damping a closed string under bounded load applied to a fixed point 
leads to nontrivial issues closely related to the above discussed ones. 
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